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Abstract 

In the context of the Miiller- Israel-Stewart second order phenomeno- 
logical theory for dissipative fluids, we analyze the effects of thermal 
conduction and viscosity in a relativistic fluid, just after its departure 
from hydrostatic equilibrium, on a time scale of the order of relax- 
ation times. Stability and causality conditions are contrasted with 
conditions for which the "effective inertial mass" vanishes. 
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1 Introduction 



The study of the evolution of self-gravitating systems (even in the spheri- 
cally symmetric case) requires the use of numerical procedures and/or the 
introduction of simplifying assumptions. The former usually lead to model 
dependent conclusions and the later are frequently too restrictive and/or 
deprived of physical justification. 

An alternative path to this question consists in perturbing the system, 
compeling it to withdraw from equilibrium state. Evaluating it after its de- 
parture from equilibrium, it is possible to study the tendency of the evolution 
of the object. This is usually done following a first order perturbative method 
which neglects cuadratic and higher terms in the perturbed quantities. If the 
relevant processes occuring in the self-gravitating object take place on time 
scales which are of the order of, or smaller than, hydrostatic time scale, then 
the quasistatic approximation fails (e.g. during the quick collapse phase pre- 
ceding neutron star formation). In this case it is necessary to evaluate the 
system immediately after its departure from equilibrium, where immediately 
means on a time scale of the order of relaxation times. 

This approach has proven to be useful in the dissipationless case Jl|, |2], ||] . 
Nevertheless, it has been recently found 0, that for non- viscous dissipative 
case this approach cannot always be applied. In particular, the goodness of 
the first order perturbation method can be examined by means of the value 
of the local parameter 

kT 
r(p + p) ' 

where k is the thermal conductivity, T is the temperature, r is the relaxation 
time for thermal signals and p and p are the energy density and the radial 
pressure respectively. If a = 1 then it has been shown that ||, |] the effec- 
tive inertial mass density of a fluid element vanishes and becomes negative if 
a > 1. The point for which the system reaches condition a = 1 is called the 
critical point. This strange behaviour of matter, at, and beyond the critical 
point might suggest that first order perturbation method fails under such 
conditions. In some cases, systems with a value of a close to, or beyond, the 
critical point are forbidden by causality conditions, but this is not always 
true. The existence of this critical point seems to take a special relevance 
for the last ones. Efectively, if causality conditions do not forbid the criti- 
cal point, then there can exist systems that cannot be studied using a first 
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order perturbation method. Furthermore, causality conditions have been 
found by means of a perturbative method up to first order in the perturbed 
quantities. Thus, for such systems, causality conditions must be taken with 
care. 

The aim of this paper is to elucidate the existence of a similar critical 
point in dissipative viscous systems. To do this we shall assume that, initially, 
the system is either static or slowly evolving along a sequence of states in 
which it is not only in hydrostatic equilibrium, but also thermally adjusted. 
Then, we shall perturb the dissipative flows and radial velocity (as seen 
by a Minkowskian observer). As it has been mentioned above, the system 
must be evaluated on a time scale of the order of relaxation times after 
the perturbation takes place. Thus, the properties of the system are still 
the same, and only time derivatives of perturbed quantities have changed 
appreciably, but not the quantities themselves. 

In order to find the critical point we shall use the transport equations for 
the dissipative flows and the radial momentum conservation equation. The 
Eckart-Landau transport equations [[?], |8| imply a vanishing relaxation time 
for dissipative flows. The adoption of these equations is not advisable here 
for two reasons: First, they predict an infinite speed for thermal and vis- 
cous signals propagation and unstable equilibrium states || . Second, we are 
evaluating the system immediately after perturbation (in the sense described 
above). Thus, to be consisten with this choice we must use transport equa- 
tions with non vanishing relaxation times. In this work, transport equations 
are introduced using the Muller-Israel- Stewart second order phenomenolog- 
ical theory for dissipative fluids ||, |10|- After the condition for the critical 
point is found we contrast stability and causality conditions with this one. 
We shall see that, in some cases, the critical point may be reached without 
violating causality conditions and, as it has been mentioned above, causality 
conditions must be used with caution. Finally, we show that neutrino trap- 
ping during gravitational collapse 0, |TJ, [TB| can lead to values of a beyond 
the critical point. 

The paper is organized as follows. In the next section the field equations, 
the conventions, and other useful formulae are introduced. In section 3 we 
briefly present transport equations. In section 4 the full system of equations 
is evaluated at the time when the object starts to depart from equilibrium 
finding the expression for the critical point. Finally, a discussion of the 
reliability of stability and causality conditions is given in the last section. 
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We adopt metrics of signature -2 and geometrised units c = G = 1 
throughout the text (except in the example preesented in last section). The 
quantities subscripted with a denote that they are evaluated at the surface 
of the sphere. 



2 Field Equations and Conventions 

We consider spherically symmetric distributions of collapsing viscous fluid, 
undergoing dissipation in the form of heat flow, and bounded by a spherical 
surface S. 

In Bondi coordinates |TJ| the line element takes the form 



ds 2 = e 2/3 



V 



du + 2dudr 



- r 2 (d9 2 + 



sin 



where u = x° is a timelike coordinate (g uu > 0), r = x l is a null coordi- 
nate (g rr = 0), and 9 = x 2 and (p = x 3 are the usual angle coordinates. 
The w-coordinate is the retarded time in the flat space-time and therefore, 
u — constant surfaces, are null cones open to the future. V and (3 are functions 



of u and r, and the "mass function", m(u, r), can be defined as [[14 

V = e w (r - 2m(w,r)). 



(2) 



Bondi and Schwarzschild coordinates (T, R, 0, <3>) are related by means 
of the expressions 

T = u+ f r L dr , (3) 
J o V 

and 

R = r, = 6, $ = y?. (4) 

On the other hand, local Minkowskian coordinates (t, x, y, z) are related to 
Bondi's radiation coordinates by 



dt 



'jdu + ^dr 



(5) 



dx 



< \lydr, 



(6) 
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dy = rd9, 



dz 



r sin 



(7) 
(8) 



For a local Minkowskian observer comoving with the fluid, the stress- 
energy tensor splits into three terms. First, an anisotropic material part, 



frnat = ^mat + p rnat^ fj^ _ ^at^ + ^mat _ p rnat^ ^ 



(9) 



where , — St, s M = 5* p mat denotes the material energy density, p mat 
refers to the material pressure for this observer, and p™ 1 = p mat + p^ at is 
the material part of the tangential pressure. Therefore, p^ at refers to the 
material anisotropy. The second one is the radiation term, which in this 



lagrangian frame reads [13, 15 



T, 



rad 



( P 





V o 



rad q 

q fpTdd 

o 








^.rad 
Pi. 





\ 




p rad J 



(10) 



where p rad denotes the radiation energy density, q the heat flow, p rad the 
radiation pressure, and p r f d = (^p rad — p rad ^j /2 . Finally, the viscous part can 



be written as 



rpviS 
fll> 



(11) 



where h^ v = i]^ — U^U V is the spatial projection tensor, and II is the bulk 
viscous pressure. The traceless viscous pressure tensor takes, for this 
observer, the form 
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( 

7T 



V o o 






(t/2) 




\ 







(*/2) J 



(12) 



where it is the shear viscous pressure. 

Thus, for a local observer comoving with the fluid the stress-energy tensor 
in local Minkowskian coordinates is 



T,u = T™ + T™ u d + Tp = (p + P ± )U,U U - P ± r,r + (P - P ± )s,s u + 2q {fl U u) , 

(13) 
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p mat + p rad ^ 



P = p + n + il (15) 



where g M = — q5*, 
is the total energy density, 

is the radial pressure, 

p=p mat +p rad ^ 

is the non-viscous radial pressure, and 

p± = pT'+pT" - ^ + n = P± -| + n (17) 

is the tangential pressure. These physical variables are obtained as measured 
by this Minkowskian observer, and the effects of gravitation are introduced 
by means of the local coordinate transformation (A£) between Minkowski 
coordinates and the Bondi ones (|||f). The dynamics of the system (the 
radial velocity of a fluid element as measured by a Minkowskian observer 
at rest in Bondi coordinates, w) can be studied applying a Lorentz boost, 
L v J—w), in the radial direction to T Mi ,. Thus, the stress-energy tensor as 
measured by an observer using Bondi coordinates, with a radial velocity, 
with respect to the matter configuration —it?, is given by the expression [f[4"| 



i.e. 



with 



T, u = L«(-w)l4(-w)AZA$f 7S , (18) 

Tfjtv = {p + P±_)U^U U - P ± g^ + {P- P±)s fl s u + 2q {fM U v) , (19) 

s, = (20) 

^-^U^ A ^ + 6^'), (21) 



_p 1 -.. / r 1 1 — w „„ / V w 



and 

^'^Vt7^77CT>- (22) 

Note that 

U% = 0. (23) 
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The traceless viscous tensor is, for this observer 



7T 



fjbv 








—w 
1+w 



7T 
7T 

-r 2 vr/2 









• (24) 



-r 2 sm 2 (6)n/2 J 
Thus, Einstein equations for the line element flU), read 



e 2/3 r 



+ - (e 2(3 -Vt + 2f3 A V 



1-w' 



-2/3 



p + 2wg + Pw 2 ) , 
(25) 



37T r* 



(e 2f3 - V A + 2P tl V) = - (p-q(l-w)- Pw) 

v J 1 + w 



2n r 2 



A: 



1 — w 
1+w 



{p-2q + P) 



(26) 



(27) 



and 



-23 



2 Aoi - ^2 K" - 2 Ai^ + 2r (/5, u y + P tX V tX )\ ) = -8vrP ± , (28) 



where subscripts , and , 1 denote partial derivative with respect to u and r 
coordinates respectively. From (|25|) and (^), it follows 



e 2/3 r 



1 - 



+ g (l + w 2 ) + Pw) . (29) 



Next, from the conservation equation = 0, we obtain after long but 
simple calculations 



e- 2 ^ + P 1 + 
2irr 



P + p 



where 



and 



1 — 2m /r 
1 



4nrP + — 



P 



1 + w 
1 

' 1 + w 



\—wp — q (1 — w) + P] 
[p — q (1 — w) — Pw] . 



-(P±-P)--(P-P) = 0, 

(30) 

(31) 



(32) 
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Finally, taking the ^-derivative of @7|) and using (@), (|IJ), (0), (g5j), 
and 
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P + p 



2r(l - w) 



(1 + (r — 2m) 



27rr y 1 — 2m /r J 
where 

V' = P + P + n + 7r — 2g 
and from (|25|), (0), and (|26"D, we obtain 

1 



1-w 2 



+ 



V - 2 m) 



(33) 
(34) 



ran 



-4nr(r — 2m)e 



2/? 



1 - 



(pw + g(l + w 2 ) +Pw) . (35) 



All these equations will be used in section 4. 

Outside fluid distribution the metric is the Vaidya one JTH], a particular 
case of the Bondi metric with (3 = 0, and V = r — 2m. The continuity of the 
first and second fundamental forms across £ leads to the well-known result 
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q a = Pa, 



or equivalently (see [18[ for details) 

Pa = -W a p a . 



(36) 
(37) 



3 Transport Equations 



As we mentioned before we shall use the Miiller-Israel-Stewart second order 
phenomenological theory for dissipative fluids [K|. Although it may be 
not reasonable in some situations, we shall assume here for simplicity, that 
there is not viscous/heat coupling {i.e. an = «i = in Jl0|]). Thus, transport 
equations read pi 



r K Kq v + q» = kW v (T„ - TU V ) - \kT 



r c n + n 



-ce - -ct 



T K U a 
T C U aS 



q» + r K u^q v , (38) 



n. 



(39) 
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and 



1 (t U° 

T v h%hP7r a p + = 2i 1 a liv - -r]T i-^- 
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2r„<^)a> ( 40 ) 



where h^ u is the projector onto the three space orthogonal to If 1 , uo^ v = 
h"h^,U[ a; p] is the vorticity, = Uft is the expansion scalar, and n, (, and r\ 
denote the thermal conductivity, and the bulk and shear viscous coefficients 
respectively Also, T, r K , r^, and r v denote temperature and relaxation times 
respectively. Overdot denotes A a p_ = U x A a p...-\ ) and the shear tensor is 
given by 

1 



a 



(41) 



The traceless viscous pressure tensor tx^ u is given by the expression (p4"|). 

Observe that, due to the symmetry of the problem, equations (|38D , and 
([40]) only have one independent component . 

Let us now write the expressions for different terms in (|38|) , they are 



-2/3 



l + W 

+r K qw 



v w 

r 1 — w 



Q.i 



r K qwe- 213 /V" 
1 — w \r 



-23 



1 - W 2 )(l + w) 



V w 

Wn H Wa 



T K qw 



1 



2m 

r 



r 1 — w 



t r V I 



q = qe 



r y/T 



nh ru T u 



KC 



-20 



l + W 



To 



V Ti 
r 1 — wJ ' 



nTh ru U v = nTe 



-20 



V 



1 — w 



1 

— I 

2 

-nT 2 q 

2 H l + w 



kT 2 



Qq r 



kT 2 



+ 



7 / w 



r \l — wj \kT 2 



(43) 

(44) 
(45) 



(46) 



9 



and 
where 



T K uTq v = r K h ra h^U [am q u = T K h ra U [a . A qP = 0, 



1 + w \2r 



1 



7T-h-7& + 77 



r /1—w 



2VJ VKl+w 



1 



(i + wy(i -w) 

and the expansion scalar is given by 



1010,1 + yiX - w)w fi 



e = e - 



/3 / r v 3/2 



r VI - w 2 Vr 



+ Ai + e 



r iw 



J 



1 — w 



For (15^), we have 



1 + w 



r Vl - w 2 (1 + io) 



m m i 



V a/1 — \ r2 r 



it) i r 

~ 77 w ,0 



1-w V 



r 1 — w 



VV 1 + w 



V w „ 

n, + -i n,i 

r 1 — w 



and 



V W { T{ 



Finally, for the different terms in ( fl0"D we get 

...2 



Tree = ^" 7r ) 



T v hghgTT al3 = T v U a ne S]C 



2 6 Tr> 



r 1 — w 



2r]a ee = 2r] 



Vy 1 + w 

r w 



V w 

7T,0 + — Z 7T,1 

r 1-w 







— e 



r VI - »' 



=r H r 

2 3 
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< / \a 



4 ' V V\ 1 + w 



+ 



V W ( T„ 



r]T J r 1 — w \T)T J 1 



(55) 



and 



2T v irf e uj e)a = 0. (56) 
Note that for this observer the shear scalar and vorticity scalar are given 



by 



and 



a 



\ 3 



r 2 - e-f** 



'V w 



r vT 



= h^h&UteflUw = 0, 



(57) 



(58) 



respectively 

Transport equations (|38] - f40D , together with (|30D will be evaluated after 
the system departs from equilibrium, neglecting terms of order 0(w 2 ) and 
higher. 



4 Departure from Hydrostatic Equilibrium 

We assume that, before perturbation, the system is slowly evolving along a 
sequence of states in which it is close to hydrostatic equilibrium and ther- 
mally adjusted - the so called complete equilibrium |21], p. 66]. Thus, the 
radial velocity, as seen by a Minkowskian observer, is small. This means that 
cuadratic and higher terms in w may be neglected in a first order pertur- 
bation theory. A system is thermally adjusted if it changes its properties 
considerabily only within a time scale T c ha that is large as compared with the 
Kelvin-Helmholtz time scale tkh- Thus, before perturbation we can assume 
that the w-derivatives of the perturbed quantities can be neglected up to first 
order, and consequently 

q t0 ~ n ~ 7r ~ w ; o ~ 0{w 2 ). (59) 

On the other hand, the hydrostatic equilibrium can be justified in terms of 
the characteristic times: If the hydrostatic time scale r hyd ~ Jr 3 /m is much 
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shorter than the Kelvin-Helmholtz time scale tkh ~ m 2 /2rl, then inertial 
terms in the equation of motion Tft = can be ignored. This condition will 
be accomplished for small values of luminosity I, and consequently for small 
values of q. Thus, before perturbation, we can assume q ~ 0(w) in the whole 
system. It seems also reasonable to assume that in such system, bulk viscous 
pressure and shear viscous pressure must be also small (i.e. II ~ tt ~ 0(w)). 

Note that from (|29f) and (|35|) /3 , Vq and m are of order w. Therefore, 
their products and second time derivatives may be neglected (0(w n ); n>2) 
(an invariant characterization of slow evolution may be found in lfffi| ). 

We shall evaluate the system immediately after perturbation (in the sense 
described in the introduction). Physically, this implies that the perturbed 
quantities (u, q, n and II) are still much less than unity. Nevertheless, the 
system is departing form hydrostatic equilibrium and thermal adjustment. 
Thus, the w-derivatives of the perturbed quantities are small but different 
from zero (i.e. g j0 ~ w ,o ~ tt,o ~ n ~ 0(w)). 

Thus, our initially slowly evolving system is characterized by: 

1. Before perturbation 

P,o ~ P,o ~ P±,o ~w~g~7r^n^m >0 
w }0 ~ q >0 ~ tt i0 w n j0 w 0(w 2 ) 

2. After perturbation 

P,o ~ P,o ~ P_l,o ^w^q 

In both cases we have 

P = P - ( W p + wp + q) + C(w 2 ) 

p = p — (wp + wp + q) + 0(w 2 ) 

The initially static case can also be considered. This system 
by: 

1. Before perturbation 

w = q = 7r = U = rn.fi = 
P,o = P,o = P±,o = w,o = 9,o = tt,o = fl o = 0. 



°^ (60) 



7i m II m i0 ~ O(w) 

n « c(w)' 



(61) 



(62) 
(63) 

is characterized 



12 



2. After perturbation 

w = q = tt = n = m o = p, = p,o = P±,o = , * 

io « g « 7r >0 ~ n 7^ (small). ^ ' 

And 

P = P, p = p. (66) 

Let us now start by evaluating ( |38|) . In the static case we obtain before 
perturbation 

"r - 2 7"2y" /? ' 1 ' (67) 

and immediately after perturbation, neglecting terms of order 0(w 2 ) and 
higher, 

r K q fi = -kTw, + O(w 2 ). (68) 

For bulk viscous pressure equation, we obtain in the initially static case after 
perturbation 

r c Tl fi = (w fi + 0(w 2 ). (69) 

Expressions (|68"D and ( |69D are also obtained for the initially slowly evolving 
case applying ( |S0|) and dOTf) - 

The evaluation of the equation for the shear viscous pressure ( fiOD yields, 
after perturbation, for both possible initial configurations 

4 

r„7r, = -vw :0 + O(w 2 ). (70) 



Thus, the three transport equations (|38| - ^0D , evaluated after perturbation, 
lead to expressions which are the same for the two initial configurations. 

Finally, let us evaluate conservation equation = ( |3Tj) after pertur- 
bation. In the initially static case, we have before perturbation condition 
). Therefore, equation Q3"0D becomes 



m 

4nrP + — 



l(p±-p) = 0, (71) 



1 — 2m/r 

which is the equation of hydrostatic equilibrium for anisotropic fluids, and 
—R denotes the total outward force acting on a given fluid element. After 
perturbation we obtain, using (|33|) , Q35|) , and ( |68H70| ) 

- R - 2 SM {1 ~ aUw ^ (72) 
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where 



or equivalently 



1 



a 



(p + p) \ 2r c 3r »? T * J 
e*R 



(1 -2m/r) j 

— 7 r — x 1 - a . 

(P + P) 



(73) 



(74) 



Assuming the second initial case (slowly evolving), conservation equation 
T£ = with (|33|) , Q35|) , and conditions (|68| - [70|) lead, before perturbation, to 
expression 



-2/3 



(r — 2m) 



(p + p) - 8vrr 2 (p + p) e w (pw + q + pw) +R — (P - P) 

(75) 



where 



R = P X 



P + p 



1 — 2m jr 



m 

4nrP + — 



(76) 



and — F (as well as —R in the previous case) may be easily interpreted as 
the total outward force acting on a given fluid element. After perturbation 
we obtain from (RBI) and (BF) 



0, 
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2nr (r — 2m) 
(p + p + Tl + n- 2q) -2 (p + p) (wfi + 4:irre 2l3 (pw + q + pw 



(77) 



Using this last expression together with 
after perturbation 



„ 2e- 2l3 (p + p) r 
F = — )!L , \l-a]xw 



[1 — 2m/r) 



we obtain for equation (B0|) 

(78) 



or 



w, 



,20 p 



[1 - 2m/ r) 



x [1 — a] 



-i 



(79) 



(p + p) 

Equations (0) and (|79|) may be compared with the Newtonian form 
Force = mass x acceleration, 

where here the term 
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1;) [!_„], (80) 

(1 - 2m/r) L J v ; 

stands for the effective inertial mass . This one vanishes for a = 1, implying 
the vanishing of —F, even though the time derivative of the radial velocity 
is different from zero. As it has been mentioned in the introduction, a = 1 
corresponds to the critical point. This one coincides with the given in [|J if 
( and r] are zero. 

Note that the effective inertial mass decreases as a grows. As a ap- 
proaches to unity the system seems to be more unstable, and for a ~ 1 
a vanishingly small radial force leads to non zero values of w t Q. This fact 
contradicts the assumption that the hydrostatic equilibrium corresponds to 
a vanishing total radial force , and consequently the reliability of a pertur- 
bative approach is in question under such condition. This approach also 
predicts an anomalous behaviour beyond the critical point. If a > 1, then 
an outward force (-F > 0) implies an inward acceleration (u> o < 0). Thus, 
we may conclude that we can neglect cuadratic and higher terms in the per- 
turbed variables only if a is not close to, or beyond, the critical point. 

Causality and stability conditions || have been found using a perturbative 
method up to first order. Therefore, it seems interesting to answer to the 
following question. Are systems with a ~ 1, or a > 1, always forbidden 
by causality conditions? If the answer is no, then the reliability of causality 
conditions is uncertain for such systems. In the next section we shall try to 
answer this question. 



5 Discussion 

In order to contrast the a ~ 1, and a > 1 conditions with stability and 
causality conditions, it is convenient to write (|73|) , using the notation adopted 
in ||. Thus, 

According to linear perturbation theory 0, causality and stability re- 
quires 
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2& + /3 

and 



(P + P)>^ r , (83) 



where 

if = 1 — (85) 
tic v Kt 

The adiabatic contribution to the speed of sound is denoted by c s , n is the 
particle number density, and c v , kt, and a p denote specific heat at con- 
stant volume, isothermal compressibility, and thermal expansion coefficient 
respectively. As usually, they are defined by 

c v = T (^) > (86) 

1 (dn\ , 
1 n\dp) T ' ' ' 

and 

« P =-:b • (88) 



n 



dT 

f ) 



As we already mentioned, if the two viscosity coefficients vanishes, we 
recover the result found in H. In this case it can be shown that the critical 
point is very close to the point where ( j8^ - ]8lf) break down ||] for small values 

of4 

Let us now consider the case where there is only bulk viscosity (re = r\ — 
0). In this case fti, oo [23|, and the critical point is overtaken if 

Po(p + p)<\, (89) 
is satisfied, whereas causality requires 

(3 (p + p)>—^. (90) 

Therefore, it appears that the critical point is forbidden by causality and 
stability requirements. 
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In the pure shear viscosity case (k = ( = 0), the critical point is overtaken 



if 



Hp+p) < - 



(91) 



whereas the most restrictive causality condition, is given (in this case) by 

«P + 1>>> jjJl^SJ- («) 

Again, the critical point is beyond the point where causality is violated. 
Let us now consider the general case («, £, 77 > 0). We may write as 

'l-c 2 \ 1 nTc^ 2 



^ + 3k + S <( " + P) 



ft 2p in Tc v - 2 ' 



(93) 



We are going to find conditions for which the point where causality is violated 
is beyond the critical point. Assuming that our system is at, or beyond, the 
critical point (a > 1), then we should demand by virtue of fl8T|) and (|93|) 



or, after some elementary algebra 



l-c 



1 nTc^ 2 
ft ~ 2ftnT Ct) - 2 : 



f3\nTc v 



2 



1 + C 2 



(94) 



(95) 



and combining it with (p|) 



^ + P / a t , \ / P + P , 



nTc„ 
implying 

which is equivalent to 



1 - 



nTc v 

1 __ 

j3inTc v 2 



1 __ 
PinTcy 2 

> 0, 



1 + c 2 



ft > 



1 



nTc„ 



2-K 2 



(96) 



(97) 



(98) 



Thus, if K 2 < 2, it is, in principle, possible to attain the critical point without 
violating causality conditions and (|84|). Note that condition (K3j) has not 
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been used in this calculus, so it should be demanded in addition to condition 
K 2 < 2. Therefore, from fl83l) , fl85l) , and (|98|) , the critical point can be 



overtaken without violating causality conditions if the following conditions 
are satisfied 

l__2e_V < 2 (99) 

TIC V K,T/ 

kT + / , < 1 (100) 



r K (p + p) r v (p + p) 



Note that conditions ( |9"9"D and ( p.00| ) do not imply necessarily a > 1, but if 



causality conditions and a > 1 are accomplished, then conditions (|99|) and 
( |100[ ) must be fulfilled. An example of this situation is an ultrarelativistic 
monoatomic ideal gas (p ~ 3nkT/2, p « n/cT, 7 = mc 2 jk^T <C 1). For this 
fluid a p = T _1 , c„ = 3/c/2, k t = p -1 , « = 4T _1 pr K /5, £ = 7 4 p r c/216 and 
i] = 2pr v /3. Then (|99"D and ( |100|) are accomplished, but a ~ 0.5. 

On the other hand if K 2 > 2, then the critical point is less restrictive 
than (|3lD, and causality and stability conditions can be used freely. Thus, 
we face the following alternatives: 

1. In the case of pure bulk or shear viscosity (without heat conduction), 
the critical point is well beyond the point where causality breaks down. 
Therefore, the system should not reach the critical point in those cases, 
and linear approximation can be applied to find stability and causality 
conditions. 

2. In the non-viscous case, the critical point is very close to the point 
where causality is violated for small values of the sound speed. Since 
the linear approximation is not reliable close to the critical point, then 
it might be possible for a given system to attain the critical point. 

3. In the general case it may happen than causality breaks down beyond 
the critical point. Thus, it appear that there exist situations where a 
given physical system may attain the critical point and even go beyond 
it. 

It is worth noticing that condition a > 1 can be accomplished in non 
very exotic systems. One of them is an interacting mixture of matter and 
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Figure 1: Temperature for which a = 1 as a function of energy density. 
Systems with a > 1, are above the line. 



neutrinos, which is a well-known scenario during the formation of a neutron 
star in a supernova explosion. In this case the heat conductivity coefficient 



is given by [24, 25 



k = ^bT\ (101) 

where r is the mean collision time, h = 7N u a/8, N v is the number in neu- 
trino flavors and a is the radiation constant. Assuming that two viscosity 
coefficients vanishes, and p p then 

kT kT 
a = — r ~ — . 102 

r K (P + p) rp 

Using usual units, the critical point is overtaken if 



T > ~ 4.29 x l y/ 4 , (103) 
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where we have adopted r ~ r K , N u = 3, T is in Kelvin and p is given in g 
cm -3 . The values of temperature, for which a = 1, are presented in figure 
1 as a function of energy density. These ones are similar to the expected 
temperature that can be reached during hot collapse in a supernova explosion 

We would like to conclude with the following comment: for degenerate 
matter, when thermal conductivity is dominated by electrons, thermal relax- 
ation time may be of the order of milliseconds (or even larger), due to larger 
mean free path of electrons [[27] , but this is of the same order of magnitude 
as the time scale of the quick phase preceding neutron star formation. There- 
fore for this last scenario (at least) , the basic assumption of our approach is 
justified. 
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